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Abstract — In this paper we present a new analysis of energy 
consumption in cellular network and we focus on the distribution 
of energy consumed by a base station. We first define the energy 
consumption model, in which the consumed energy is divided 
into two parts: the additive part and the broadcast part. The 
additive part is served to communicate with different users while 
the broadcast part is served to transmit the same information 
to all users in the cell. We then model user's activity as ON- 
OFF process in time. We also model user's mobility as a random 
process and we define the high mobility regime in which users 
move very fast. We are able to provide analytical expressions 
for statistics of consumed energy and bounds on the distribution 
of energy in motionless case. We then consider the impact of 
mobility and we mobility reduces moments and central moments 
of the additive part. We also show that high mobility regime, the 
variance of energy tends to zero. This is a strong result as it holds 
true for almost any mobility model. We are able to characterize 
the convergence decay rate of the variance in function of user's 
speed. We propose two applications of the model. The first one 
is to dimension cell radius in the economical point of view and 
the second one is to dimension cell battery in the sites that do 
not have access to electricity. 



I. Introduction 



A. Motivation 



Cellular communications have had a phenomenal progres- 
sion due to recent technological advances in both cellular 
networks and cellular telephone manufacturing. They will 
experience even more growth in the next decade. Data traffic 
is becoming more and more dominant in mobile networks. For 
the first time in history, the volume of worldwide mobile data 
traffic exceeded that of voice traffic in December 2009 (HI). 

Such demands, on both number of subscribers and traffic 
capacity, are driving forces for the telecommunication in- 
dustries to develop new technologies, and for operators to 
invest more and more in infrastructures. According to several 
studies, 0.5% of world-wide electrical energy is responsible 
to cellular network and 80% of this is consumed at base 
station sites. Up to 90% of cellular energy consumption is 
operator's operational expenditures (OPEX) (|2],f3),|D). A BS 
connected to electtical grid may cost approximately 3000$ 
per year to operate, while an off-grid BS may cost ten times 
more. Considering the rising energy consumption of mobile 
networks, it becomes clear that energy cost is critical for 
operators' OPEX. The expense for energy to run cellular 
networks is expected to triple over the next seven years. 



Consequently, one of the main objectives for new generation 
of cellular network is to reduce energy consumption and 
carbon emission and to improve energy utilization efficiency. 
That leads to the concept of "green cellular network". The 
interest of green cellular network is not limited to ecological 
reasons but also includes on economical benefits. We refer to 
1 3 1, |5| and reference therein for green techniques for cellular 
network. 



B. Contribution summary 

In this paper, we investigate different aspects, which sur- 
prisingly have not been investigated in the literature. We are 
interested in the disttibution of energy consumed by a base 
station during its operating period. As mobility is the main 
issue in cellular network, we also take it into account. More 
precisely, we model users' initial position as a spatial Poisson 
point process, users' activity as ON-OFF random process and 
users' motion as random processes. We define the energy 
consumption model in which the consumed energy consists 
of the additive part to communicate different information to 
different users and the broadcast part to transmit the same 
information to all users in the cell. We are able to provide 
analytical expressions for the moments of the additive part. 
We provide bounds on the disttibution of the consumed energy 
and we show that mobility reduces central moments and 
moments of the additive part. We also show that in high 
mobility regime, the variance of consumed energy reduces to 
zero and we characterize the convergence rate in the function 
of users' speed. We apply the achieved results to solve two 
dimensioning problems. The first one is to dimension optimal 
cell radius and the second one is to dimension base station's 
battery. 

The rest of the present paper is organized in the following 
way. In the section [TT] we describes the model, including the 
ON-OFF model and the mobility model. In the section [III] 
we provide analytical results for the case of motionless users. 
In the section [IV] we analyze the case where users move. In 
the section [V] we refine some results in the section [IV] to 
the case where users' speed is constant. In the section [VI] 
we consider two applications of the model: dimensioning cell 
radius and dimensioning base station's battery. Here is the list 
of mathematical notations appeared in this paper: 
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Symbols 


Definition 


K+ 


[0,oo) 


D 


differential operator 


P(A) 


probability of event A 


E[X] 


expectation of random variable X 


V [X\ 


variance of X 


m„ [X] 


n th on } er moment of X 


c„ [X] 


n th or( j er cen tral moment of X 


B(x,r) 


ball of radius r centered at x (in d dimension) 


1{A} ( x ) 


indicator function 


C [X, Y] 


covariance of two random variable X, Y 


Q(u) 


-J= f n e"4 da; 

\/2tt J— 00 


Q(u) 


-3= f°° e -4 dx 


g(u) 


(1 + u)ln(l + u) - u 



Table I 

Mathematical notations. 




Uplink and downlink 



Figure 1. Power consumption model. 



II. System model 
A. Model for energy consumption 

We suppose that there is a cellular network with multiple 
base stations on R d and there is a base station located at 
the origin o administering a geographical region C around 
o. Practically speaking, base stations are deployed on a plane. 
Hence d — 2. However for the sake of generality, we compute 
the formulas for any value of d. Note that considering d = 3 
can be interesting if we study the deployment of wireless 
access points in a building. We assume that there exists 
< Ri < R such that B(o,Ri) c C C B(o,R) and C is 
convex and compact. We define Rmf — ini^{C C B(o, R)} 
For a given spatial configuration of active users on M. d at an 
instance, denoted by rj, users located inside C will be served 
by o, users outside this region will be served by another base 
station (or in outage regime). The power consumed by the 
battery of the base station o can be divided into two parts: 



• the power dedicated to transmit, receive, decode and 
encode the signal of any active user. The cumulated 
power over the whole configuration is then of the form 
T^xe-q < ^( a; )' where ^ is a function to be defined later. 

• the power dedicated to broadcast messages. In order to 
guarantee that all active users receive these messages, the 
power must be such that the farthest user in the cell is 
within the reception range (if the system performs power 
control) or all the cell is within reception range (if the 
system does not performs power control). Thus, the power 
is a function of max l£ , n c \x\ where \x\ is the Euclidean 
norm of x (the power is equal to if r\ n C = 0). This 
function is constant if power control is not performed. 

It follows that the total consumed power is given by: 



P{ri) = P A (rj) + P B (r,), 



(1) 



where P A (r}) = ^xe,,^) and p b(v) = VKNI). IMI = 
max iei)nc \x\ if 77 n C ^ and || 77H = Oif?7nC = (the 
subscript A stands for "additive" and B stands for "broadcast"). 
For a very simple propagation model (without fading and 
shadowing), the Shannon's formula states that for a receiver 
located at x, the transmission rate is given by 

^io g2 (i + p e ;(x)), 

where W is the bandwidth, P e is the transmitted power and 
l(x) is the pathloss function. Generally, the function I : M. d — » 
[0, 00] takes the form l(\x\) where I : [0, 00] — > [0, 00] is a non 
decreasing function. This implies that in order to guarantee a 
minimum rate at position x, P e must be proportional to l(x). 
Thus, it is sensible to choose <fi as 

<t>[x) = a.~l(\x\)l {xeC} 

with a > 0. The function tp is chosen as 

1 / \ _ \ b.l(\x\)l{ x& c}i if power control is performed; 

1 bd(Rinf), if power control is not performed. 

We can divide models for path loss into two categories: 

• singular path loss model l(x) = K \x\^ 1 where 7 is the 
exponent path loss parameter and K is a positive constant. 

• non-singular path loss models like l(x) = K(ro V \x\) 1 
orl{x) = K{l + \x\~ 1 )- 1 . 

More generally, we make the following assumption: 
Assumption 1. The transmitted power depends only on the 



distance to the base station (j>(x) 



M)i { 



zee}- 



Further- 



more, 



and tp are continuous non decreasing function on 



We denote ip(x) = ip(\x\)li xe c}- This implies that <\> and 
ip are always bounded function. Apart from the above model 
for power consumption, we can define a power consumption 
as a general functional depending on the configuration of users 
P G :tl Rd —►[Coo]. 

If lu = {u>t,t > 0) is a process of time varying configura- 
tions, the total consumed energy between time and time T 
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is given by 

J T := J T (u>,T) = [ P(uj s )ds- 
Jo 

As previously, we also define J a and Jb by: 

Ja:=Ja(^,T)= f P A {uJs)ds 
Jo 

and 

J B := J B (w,T) = / P B (u s )ds- 
Jo 

The same definition for Jg(uj,T) if the system applies power 
consumption 

Also we denote C(r) = Cn B(o,r) and C(r) = Cn 
B(o, r). For a configuration i^, we denote Xj, the point of v 
such that \x v \ = \\v\\ (if there are more than one point then 
Xnu is randomly chosen among these points). 

B. Model for mobility of users 

In this section, we introduce the mobility models for users. 

Consider the functional space Z3(R + ,R d ) of rcll functions 
on W l equipped with the Skorohod topology(see, for example 
(6), page 369). It is well known that T>(R+,R d ) is a Polish 
space. The subset T> (R+,R d ) = {/ G D(R+, R d ), /(0) = 
0} equipped with the Skorohod topology is also a Pol- 
ish space. We consider a probability distribution Pm of 
a random variable M — (M(t),t G R + ) defined on 
the associated Borelian er— field of the space £>o(R + , R d )- 
Each realization of M can be represented as a rcll tra- 
jectory of on R d . Also, this probability is completely de- 
termined by the distributions of finite marginal distributions 
P(M(ii) G M(t n ) G .)(ti,...,t„ > 0). In some situa- 
tion, for convenience we can assume that M(t) = o for t < 0. 
M is said to satisfy the property T if P(M (t x ) = M(i 2 )) = 
for any < t\ < t 2 . 

If mobility is considered, then each user is associated with 
a mobility process on R d . We make the following assumption: 

Assumption 2. 

• Motion trajectories of users are i.i.d mutually indepen- 
dent and have the same distribution as that of M. 

• Motion trajectories of users do not depend on the initial 
position of users. 

More precisely, consider a user i initially arriving at Xi is 
associated with an independent version of M, namely Mj and 
an arrival time Tj. This user will move during its sojourn along 
Mi, i.e the position of this user at time t > Tj is x + Mi(t — 
Tj). The random process (Mj)jgN are mutually independent. 
Examples for mobility model are as follows: 

. Motionless users: M(t) = o,Wt G R. 

• Brownian motion users: M(t) = c(t)E>d(t) where c(t) G 
M is a continuous function in R + and Bd is a standard 
Brownian motion on R d . 

• Completely aimless users: M(t) — tv where the speed 
of user v is random whose direction is chosen randomly 



and uniformly over the d-dimensional unit sphere and \v\ 
is a positive random variable. 

• Combination of two above models: M(t) — tv + 
c(t)B d (t). 

• High mobility regime: let e > be a small parameter, the 
high mobility regime consists of considering the mobility 
process (M/e)(t) = M(t)/e and we want to study the 
behavior of the system when e — > 00. The high mobility 
regime of the completely aimless mobility model with 
constant speed 1 1? j is the same as considering |«| — > 00, 
i.e when the user's speed is very high. 

C. ON-OFF model for users' activity 

An ON-OFF process on the real line alternates between 
values 1 (for on state) and (for OFF state). ON-periods (and 
OFF-periods) are i.i.d positive random variables. Furthermore, 
the sequences of ON-periods and OFF-periods are indepen- 
dent. Each realization of an ON-OFF source is a rcll function. 
An ON-OFF process is called exponential if ON-periods and 
OFF-periods are exponential distributed. 

More precisely, we consider an ON-OFF sources (I(t),t G 
R) such that the ON-periods are continuous positive random 
variable of mean [i^ 1 > and the OFF-periods are continuous 
positive random variable of mean fi^ 1 > and denote by U 
and V the generic ON-period and OFF-period. We can write: 

00 

I ( t ) = ^2 1 {T2i<t<T 2l + 1 }- 
i= — 00 

where ... < T 3 < T 2 < T_i < T < T x < T 2 < T 3 < ... such 
that (T2i — T2i~i)°Z_ 00 are i.i.d and have the same distribution 
as V, and (T^+i — T2i)°^_ oa are i.i.d and have the same 
distribution as U. We can assume that T V T x > 0. I can 
be seen as a random variable on Z?(R, R) with probability 
measure dP/. We assume that / is stationary. From [7| for 
example, we have P(I(t) = 1) = -^f-^ = n 1 and P(/(t) = 
0) = — *r — = 7i"o for all t. Furthermore, we have 

roo 

P(i(t) = me[o,u)) = ^j- j P(u> s )ds. 

for all u > 0. 

We make the following assumptions: 

Assumption 3. The positions of users at t = follow a 
Poisson point process N = {Xi}i>i of intensity measure 
Ada;. User i is associated with an ON-OFF process of activity 
(Ii(t), t G M), i.e users are active during their ON-periods and 
are inactive during their OFF-periods. The activity processes 
of users are assumed to be i.i.d and have the same distribution 
as that of {I(t),t > 0). 

Following the above assumptions, the configuration of active 
users at time t is 

W * M = ^2 1 {h(t)=l}5x i +M % (t)- 

i>l 
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The system can be described by a Poisson point process on 

I d xD(R,l)xD(K+,R <J ) 



4> 



I,M 



of intensity A dx x dP/ x dPj\/. The consumed energy is 
defined in the same way as in the previous subsection for the 
time varying configuration process lu m = (u) M (t),t > 0). 
In particular, the additive part of consumed energy can be 
rewritten as: 

i-T 

iit^Xi + Mi(t)) dt 



J A (w M ,T)=Y / ( h 
i>i Jo 

and the broadcast part is: 

Mu m ,t)= / T ?(|| w f ||)dt- 
Jo 

The total consumed energy is J t (lo m ,T) = J a (lo m ,T) + 

When users are motionless, i.e Mj(t) = o, the system is 
described as a Poisson point process 

& = {(X l ,I l )} t > v 

of intensity measure A dx x dP/. In this case, we drop the 
superscript. Thus, the configuration of users at time t is: 

n 

w t — l{j ( (t)=o}£x 4 - 

z=l 

The additive part of consumed energy is 

j a (u,t)= y2<p(Xi) [ h{t) dt 

i>l J ° 

and the broadcast part is: 



Jb(u,T) 



^(|k M ||) dt- 



The total consumed energy in this case is Jt(w,T) = 
J A (u,T) + Jb(<jj,T). In the next two sections we present 
analytical results on the motionless case and the general case. 

III. Motionless case 

In this section, we assume that users are motion- 
less. We derive analytical expressions for the moments of 

J A (u,T),J B (u,T),J T (u,T) in this case. Let G(h) = 
J Rd h{x) dx for h : R d -> K. 

Theorem 1. The moments of Ja{u,T) and the central mo- 
ments of J a (uj,T) are given by: 

m„ [J A (u,T)] = 

B n (Arm [A{T)\ G{4>), Am„ [A(T)} G(4> n )) , 



and 



c„ [J A (w,T)] = 

B n (0, Am 2 [A{T)\ G(<j> 2 ), Am„ [A(T)} G{4> n )) 



In particular, the expectation of J a (oj,T) is given as: 

E[J a (uj,T)] = Avri / (j>(x) dx. (2) 

Proof: The theorem is derived from theorem[T3] see more 
in details in [8| (section 7.3, theorem 31). ■ 

We see that, from the above theorem, the expectation of 
J A (cj, T) depends on the distribution of ON -periods and OFF- 
periods only by the activity rate pof/ii. 

Applying theorem [14] an error bound for Gaussian approx- 
imation of J A (ui,T) is found as follows: 

J A (u,T)-E[J A (u,T)] 
V[Ja(w,T)] 



Theorem 2. Let J A (to,T) 
u we have: 



IP (Ja(u,T) > u) ~ Q(u)\ < 



then for any 

m 3 [A(T)} G(f ) 
V\(m 2 [A(T)]G(0 2 ))5 



The above bound decays as O as A — > oo. As shown 

in Lemma |20] T 2 > m 2 [A(T)\ > tt 2 T 2 and m 3 [A{T)\ < 
T 3 , so the bound decays as 0(1) as T — > oo. We obtain a less 
sharp bound but depending only on the activity rate ni but not 
on the distribution of ON-periods and OFF-periods and on T: 



\P (Ja(u},T) > u) - Q(u)\ < 



J Rd c/) 3 (x) dx 



As already noted in the case of exponential ON-OFF source, 
when T goes to infinity, m 2 [A(T)] ~ vr 2 T 2 and m 3 [A(T)\ ~ 
irfT 3 . Consequently, in this case the bound has the following 
limit when T — > oo: 

j Rd 4?{x) dx 



'7T1 



a (Jr<* ^ 2 ( x ) dx ) 



Theorem 3. The joint distribution of 
given by: 



'(II 



il,...,M„) = X!^ 1 ) 



E 



l<fci<...<fci<n 



,l(tfc! ,---,t k .)A(C(max{u kl ,...,u ki }))_ 



In particular, 



(3) 



The n order moment of Jb(w,T) are given by: 



m„ [J B (u,T)] 
dFn\.. II 



dt%... dt n 



[0,T)« 



ip(u 1 )...ip(u n ) 



Proof: See (section 7.3, theorem 33). ■ 
We have obtained closed form formulas for V[Jb(cj,T)] 
and V [Jt(w, T)], however it requires to solve triple integrals. 
We now present bounds for V [J B (u),T)] and V [J T {u),T)\, 
which requires only one easily computed integral. 



Theorem 4. We have: 

V[J B (u,T)] < Am 2 [A(T)]G(^ 2 ), 



(4) 
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and 

Am 2 [A(T)\ G{4> 2 ) < V [J T (u, T)] < Am 2 [A(T)\ G{{4> - 



Mobility reduces moments of J A ; is 



m„ [J A (^ M ,T)] <m n [J A (o;,T)] 



Proof: The theorem is derived from theorem [16] and 
theorem [TTl see [8| (section 7.3, theorem 35). ■ 
The following theorem gives an upper bound on the distri- 
bution of Jt(uj,T): 



Theorem 5. Assume that 
let 



x) + Tp(x) < K for all x G 
a 2 = m 2 [A(T)} [ (i/j(x) + 4>{x)) 2 dx 



then 



P(J r (w,T) > E[J T (w,T)] +«) < 



exp 



A 2 c 



-.9 



/ uTK\ 



for all u > 0. 

Proof: The theorem is derived from theorem [18] see also 
[8| (section 7.3, theorem 36). ■ 

IV. Impact of mobility 

We have, by the lemma [T9j when users move, the spatial 
distribution of users remain the same as without mobility. 

Following the previous lemma, we can prove that in both 
cases with and without mobility, the consumed energy has the 
same expectation, which is somewhat surprising. 

Theorem 6. For any power allocation policy Pq, and for any 
mobility model M, the expectation of energy consumed is the 
same as in motionless case, i.e: 

E[J G (^ M ,T)] =E[J G (w,T)]. 

In particular V[J a (lo m ,T)] = E[J A (u),T)}, 
E[j B (w M ,T)] = E[J B (u,T)} and E[J t (lo m ,T)] = 
E[J T (w,T)]. 

Proof: See J8) (section 7.4, theorem 38). ■ 
The distribution of consumed energy in the mobility case is 
clearly not the same as in the motionless case. Let 



E 



f{x + M(t))I{t) dt 



In the following theorem, we characterize the impact of 
mobility on the distribution of the additive part of consumed 
energy: 

Theorem 7. The moments of J A (u M ,T) are given by 

m„ [J a (lu m ,T)] = 

B n (XFf(T, 1), \Ff(T, 2), \Ff(T, n)) 



and 



c„ [J a {lu m ,T)] =B n (0,XFf (T,2),...,XFf (T,n)) 



and 



c„ [J a (uj m ,T)} <c n [J A {u,T)]. 
Furthermore, J A (u> M ,T) is Laplace-smaller than J a (uj,T), 



E[exp{aJ4w M ,T)}] < E [exp {aJ A {w, T)}] (5) 

for all a € K. If M has the property T then the central 
moments of J A goes to in high mobility regime; i.e 



J a (uj m /\T) 
Ja(u M/ \T) 



-> 0, 



7TlA 



dx 



(E[J A ( Wj T)]) n 



as 0-. 



Proof: The expressions for the moments and central mo- 
ments come from theoremQj] The Laplace-smaller property is 
proved by the fact that e au is a convex function for all a£l. 
The moments reducing property and the convergence to is 
derived thanks to lemma [22] For more details, see [8| (section 
7.3, theorem 41). ■ 
The following theorem gives impact of high mobility on the 
distribution of J B {u M , T) and J t (lo m ,T): 

Theorem 8. The variance of Jb{u m ,T) and Jt(w m ,T) are 
bounded as follows: 

AT)] 



V[J B (. 



< 



(T,2) 



\f¥(T,2)<V[J t (lo m ,T)} < XFf +i) (T,2). 

Therefore, if M has the property T then in high mobility 
regime, the variance of Jb and Jt tends to 0, i.e 



J B (u M, \T) ,V J t {u u '%T) 



M/e 



as e -» 0- 



Proof: The theorem is proved thanks to theorem [17] and 
lemma [22] See also [8| (section 7.4, theorem 42). ■ 

The above results say that, when users move the total 
consumed energy by a base station does not change in average, 
and the moments and central moments of the additive part are 
reduced. Moreover, when users move very fast, the consumed 
energy during a time period is almost constant. We can see 
this fact as a consequence of weak central limit theorem. 
When users move faster, the configuration of users takes more 
"value" on fl R during a same period of time, thus converge 
faster to the mean. 

We find an error bound for Gaussian approximation of 
J a {uj m ,T) as follows: 



Theorem 9. Let J a {lo m ,T) = V [j a (^.t)} 



then: 



P(j A (u: M ,T) >u) -Q(u)\ < 



Ff(T,2) 



VA(V^(T,3) 



3/2 
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Proof: This is a consequence of theorem [T4l See also |8| 
(section 7.4, theorem 43). ■ 

Theorem 10. Assume that (f>(x) + ip(x) < K for all x g R d , 
then 



P (J t (w m ,T) > E [J t (lu m ,T)] + u) < 



T 2 R 2 



uTK 



eXP \ *Ff + 4T,2) 9 {f¥^(T,2) ) j 

for all u > 0. 

Proof: This is a consequence of theorem [T8l See also |8| 
(section 7.4, theorem 43). ■ 

V. Completely aimless mobility model 

In this section, we consider the completely aimless mobility 
model M with constant speed, i.e M(t) = tv where \v\ = 
constant and the direction of v is uniformly distributed. The 
following theorem characterizes explicitly the convergent rate 
to of central moments of J a, as well as the variance of Jb 
and Jt in function of user speed \v\. 

Theorem 11. Consider the completely aimless mobility model 
M(t) = vt with constant speed \v\ then 

1) V [J a (uj m ,T)] andV[j T (uj M ,T)] decay as Q(±) as 
\v\ — » oo and decay as 0(T) as T — > oo . The coefficient 
of variations CV{J a {uj m ,T)) and CV(J t (uj m ,T)) de- 
cay as 0(— -7=) as \v\ — > 00 and 0(-i=) as T->co. 

2) The error bound o£ Gaussian approximation in Theorem 
^decays as 0(|v| ) as \v\ — > 00 and decays as 0(-3=) 
as T — > 00. 

Proof: From the assumption Q] we have <fi(x) < ci for all 
x G C, 4>(x) > c 2 for all x £ C/B(o, ^-) and < ip(x) < c 3 
for all x £ C for some finite positive constants c%, o%, C3. Now, 
according to Theorem [HJ 

\F¥(T,n)=V[j A (u M ,T)] 

<V[J T ( W M ,T)] < XFf +il (T,n). 

The results then follows from Lemma |231 ■ 
We also see that, if \v\ is small, the variance of Jt(oj m ,T) 
is proportional to ^ while in the motionless case, it is 
0(T 2 ). We also notice that mobility makes the Gaussian 
approximation of J a more accurate when T is large as the 
bound decays as 8(-j=) instead of 0(1) in the motionless 
case. In the motionless case the position of users are always 
fixed over time, only their state change, thus the configuration 
of active users can take only some possible values on f2 K . 
So it is intuitive that one cannot guarantee that the Gaussian 
approximation is good if T is large. On the contrary, in the 
mobility case, the configuration of active users can take all 
possible values on f2 K . When T grows larger, it take more 
values. Thus, Gaussian approximation is better when T grows 
larger. Quite surprisingly, when \v\ is large, the variance of 
J A tends to but the bound on Gaussian approximation does 
not decrease. 



Remark 1. From the above proof and using the properties 
of Bell polynomial, we can prove that c 3 \ J a (uj M ,T)] 



e(^), Ci [j A (u 



M 



T) 



V \v\ 



-), c 5 [J A (u; M ,T)] 



VI. Applications 



A. Dimensioning optimal cell size 

Consider an operator aiming to design the optimal cell 
radius R to cover a region of total area (volume) S € M. d . 
We assume that the cells are circular. The average total cost 
of the network is assumed to be the sum of the operation 
cost during the life time of the network (says T) and the 
cost of facilities (base stations). We assume a fixed cost of 
base station regardless of its transmission range. The number 
of base station is then proportional to -Jj-, say -^j, so the 
installation cost of base stations is with ci > 0. The 
operation cost is assumed to be proportional to the consumed 
energy. 

We assume that l(x) — K. |x| 7 , <fi(x) — a.l(x)ls xe c}> 
ip(x) — b.l(R) with K > and furthermore 7 > d. From 
results presented in previous sections and from lemma [24] 
the expectations of the two parties of energy consumed by 
a singe base station are ~E[Ja(ui,T)] = XTa K y+d and 
E [J B (w, T)} = Tb'ffl with a, b' > 0. The energy consumed 
by the network during its operating time is: 



C 



7— d 



-jp {a'\TR d+ ~< + b'TR~<) = a x \TR? + b\TR 

This is a increasing function of R, which means that small 
cell systems will consume less energy than larger cell system. 
The average total cost for the network is then 



Cost(R) = axXTW + biTR 1 ~ d + 



R d 



It is interesting for operator to find the optimal R in order to 
minimize the cost function. As Cost(R) > for all R > 
and lirrii^o Cost(R) = linii^oo Cost(R) = 00 there exists 
a minimum for Cost. By differentiation, the optimal cell radius 
is the unique positive solution of the following equation: 



ai\jTR2p t d 



•&iT(7- 



As the RHS is increasing in T, the optimal cell radius 
R op t must be a decreasing function of T. This reveals a 
characteristic of the optimal choice of cell radius. In the 
economical point of view, to operate a network with longer 
life time it is preferable to exploit smaller cells system. 

If bi — 0, i.e the broadcast part of transmitted power is small 
comparing to the additive part, then the problem reduces to 
minimizing aiXTR 1 + jp-. Simple manipulations yields: 



R, 



opt 



\-fXaxT J 



That is to say theoretically the optimal cell radius is propor- 
tional to (AT)~^T3. 
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B. Dimensioning cell battery 

The proposed model can be used to dimension sites that do 
not have access to power supply facilities. In this situation, 
operators have to replace or reload base station's battery for 
each period T. We want to determine the energy level a of 
battery so that the probability of running out of energy before 
replacement(or reloading) be smaller than some threefold e <C 
1. We use results derived in the previous sections to find a. 

Assume the mobility model M and that we are able to 
estimate f¥ (T, 2) and F¥ + ^{T,n). The problem is to find 



a such that: 



P{j A (u M ,T) >a) <e. 



From theorem [10] we deduce that the following value is a 
solution: 



a\ = Tiri 



TK 



where u\ > is the unique solution of the following equation: 



K«i) 



XF¥ + JT,2)]ne 



T 2 R 2 



Now applying theorem [9J we can find another solution. 
Assume that ip(x) < Kb for all x £ C and 

F¥(T,2) 



3/2 



< e 



then J T (w M ,T) < J a (oj m ,T) + K B T and: 
P (Jt(uj m ,T) > a) < P (J a (lu m ,T) >a- K B T) 
Therefore, the following value is another solution: 



a 2 = AG(0) + u 2 y XF ¥ (T, 2) + K B T 
where u 2 is the unique positive number such that: 

F¥(T,2) 



Q{u 2 ) 



Vx(f¥(t,3) 

VII. Conclusion 



3/2 " 



Throughout this paper, we have presented a new analysis 
of the energy consumed by a single base station in a cellular 
network, we have assumed that each user is associated with 
an ON-OFF process of activity. We have derived analytical 
expressions for the distribution of energy consumed by a base 
station. We have found that, with or without mobility, the base 
station is expected to consume the same amount of energy in 
average. We have proved that mobility reduced moments of 
the additive part of energy. We have also proved that high 
mobility leads the variance of energy to 0. These results are 
strong since they hold for any mobility model. In the case 
of completely aimless mobility model, we have characterized 
the convergence rate to of the variance, which is The 
mathematical framework presented in this paper is new and 
can be served to further studies. 
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Appendix 

Poisson point process 

Let E be a cr— compact metric space (i.e E can be parti- 
tioned into a countable union of compact metric spaces) with 
a diffuse Radon measure v. The space of configurations of E 
is the set of locally finite simple point measures 

Q E := = e z< (at most countable), Zi € e\ , 

where e z denotes the Dirac measure at z € E, i.e 

5 z (A) = MzeA},AeB(E)- 

Here, simple measure means that w({z}) < 1 and locally finite 
measure means that w{K) < oo for all compact K C E. The 
configuration space Vl E is endowed with the vague topology 
and its associated er-algebra denoted by T E . 

For convenience, it is quite often to identify an element 
w = J_ e Zi with its corresponding support, i.e the unordered 
set {zi, z n }, n € N U {+oo}. Also, w(A) counts the 
number of points in A E B(E). The distribution of point 
process w is characterized by the family of finite dimensional 
distributions (w(Ai), w(A n )) where Ai,...,A n are mutu- 
ally disjoint compact subsets of E. 



8 



Definition 1. w is a Poisson point process (PPP) of intensity v 
if for all set [A\, A n ) of mutually disjoint compact subsets 
ofE: 



is nf — J E P{z) dv(z) (i = l..n). The moments and central 
moments of F are given as: 



PO(Ai) = ki,...,w(A n ) = k n ) = Yl 



MA,) WW 



ki 



kA 



mi [F] 



B, 



f(z) dv(z), f f(z) du(z),..., f f(z) du{z) 

JE JE 



If E = R d ; B is Borel algebra of R d and v{dz) = Adz, 
we will call w the homogenous Poisson point process with 
intensity parameter X on M. d . 

Roughly speaking, the number of points of w falling into 
a subset A follows Poisson distribution of parameter v{A), 
and the number of points falling into 2 disjoint subsets are 
independent. 

Linear functional 

We call F a linear functional of w if there exists / : E — > R 
such that 

f= f f( Z )dw(z) = j2m- 

We assume that / e L 1 ^). In this section we are interested 
in the distribution of F. 

Let C w (.) be the Laplace functional of w, i.e: 



(9) 



f\z) du(z),..., / f(z) du(z) ) (10) 

E JE 



C w (u) = E 



f E u(z) dw(z) 



(6) 



Theorem 12. f/9 J) The Laplace functional of w satisfies: 

C w {u) = e -/ J .(i-«-" w )d"W, utL 1 ^)- (7) 

From the above theorem the moment generating function 
(MGF) of F is expressed as follows: 



E 



S E f Aw 



gJ E 



/ B (e«*>-1) dr/( Z ). 



6 JE 



(8) 



The complete Bell polynomials B n (ai, a n ) are defined 
as follows: 



I LXJ (JU 

«p E> -E 



B n (ai,a 2 , ...,a„) 



for all a\, a n and 8 such that all above terms are correctly 
defined. 

The first four Bell complete polynomials are given as: 



Bi(oi) = ai 
B 2 (ai,a 2 ) 

B 3 (ai,a 2 ,a 3 ) = a? + 3aia 2 + a 3 

-84(01,02,03,04) = Oj + 4a^a 2 + 4aia3 + + a 4 



The Bell polynomials helps us to express central moments and 
moments of F in a simple fashion: 

Theorem 13 (Generalization of Campbell's formulas). As- 
sume that f £ n" =1 L*(-E, v). The cumulants of F = f E fdw 



and 

c t [F] = B, 

for i = 1,2, ...,n. 

Proof: ■ 
As a direct consequence, one can easily obtain from the 
above theorem two useful formulas (Campbell): 

Corollary 1. Let F = j E f dw then 

E[F]= [ f(z) du(z), feL\ v ) 

JE 

V[F]=E[(F-E[F]) 2 ] = f f(z)dv(z), f € L 2 (v)- 

JE 

We have now expressions for moments and central moments 
of F. We note that the central moments of F is always non 
negative as / is supposed to be non negative. One can ask if 
there is mean to compute the tail distribution of F. Gaussian 
approximation seems to be a first answer one can think of 
due to the central limit theorem (CLT). An error bound of 
Gaussian approximation for sum of n i.i.d random variables 
is known as BerryEsseen theorem. It is possible to find an 
alternative version of this error bound for linear functional 
of PPP. Let Q(a) = -4= P e^ 1 '/ 2 du be the CDF of a 
standard Gaussian random variable and Q be its CCDF 

Theorem 14. (gWj, 1$) Consider F = J E f dw with f E 
L 2 (y). Let F = F ~^ ] , then: 



\P(F < 0) - Q(a)\ < 



f E \f(z)\ 3 dv(z) 
(j E P(z)dv(z))y* 



(ID 



We have presented the Gaussian approximation and Edge- 
worth expansions for linear functionals. We are now interested 
in upper bounds on the distribution of F, which can be called 
concentration inequality. 

Theorem 15. (fn}, Ml/J Let M, a > 0. 

Assume that \ f(z)\ < M v—a.s and f 6 L 2 (E,v), then 



P(F > E [F] + a) < exp 



r2 - y 



a.M 



M 2 * \ V[F] 
General functional 



(12) 



We now consider a general class of functional L 2 (fl E ,P) 
(not necessarily linear). Define the difference operator D as 
follows: 

D z F(w) = F(w + e z ) - F(w), dPxdi/fl.e, 
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for F E L (fi , P), w e H fi and z e F. Note that if F = Lemma 20. We have, for all T 

f„ f(z) dz then D z F(w) = f(z). _„ . ,,„,., m „ 

Using the difference operator F/, we can bound the variance 1 — n [ \ n — 

of a functional F as follows: Proof: Since A(T) < T a.s. we have m„ [A(T)] < 

Theorem 16. flZB VF € £ 2 (fi s ,P) we have: T™. Now m„ [-A(T)] > (E [A(T)]) n = <T«by 

Cauchy— Schwarz inequality. ■ 

V \F] < E /l_D FI 2 d (13) ^ or a non ne 8 at i ve function / E and ti, ...,i n € K 

[Je Y and f(x) = if x E R d /C. Define 

Equality occurs if and only if F is linear. -*m,, \ f -^rp, ,,/, ,> „, , r/ NN1 , 

* 7 /7 *n{f,h,...,t n )= E[f(x + M{t 1 ))...f(x + M(t n ))]dx. 

Corollary 2. VF e L 2 (Q, E , P) sue/* fnaf |£> Z F| < /(z) /or 

.some non negative measurable function f : E -» K /or aW Lemma 21. We nave $* f (/, i) = J Rd /(x) di ant/: 

z, io f/zen: . 

V[F] < / f{z) dv{z)- (14) $M{f,t u ...,t n )< f n {x)dx- 

Je Rd 

Theorem 17. (i§fc£2fr Consider two functionals Moreover, if M has the property T and n > 2 fnen 

Fi,F 2 E L 2 (fi s ,P) (not necessarily linear). Assume 0» £ (/:*l: ->*n) a* e ->■ with n > 2. 

that D Z F 1 ,D Z F 2 > 0, P x v a.s. then Pmo f. See Jg] ( sect i on 7.4, lemma 39). ■ 

C [Fi,F 2 ] > 0. Lemma 22. We have: 

As a consequence, V [Fi + F 2 ] > V [Fi] + V [F 2 ]. T) < m n [A(T)] / / n (x) dx- 

Recent results concerning Edgeworth's expansion for func- 
tionals depending on Poisson point processes can be found in If M has the property T then F^ e {T, n) ->■ as e /or 
|Q3l and 0. n > 2. 

Theorem 18. /77/, Lef F E domD. Assume that \DF\ < Proof: See (S) (section 7.4, lemma 40). ■ 

M, P x v a.s., for some M > ana 1 J^, |D Z F| dv(z) < L emma 23. Let f(x) be a positive measurable function on W 1 

a , P-a.s. Then for all u > we have ^ ^ f{x) = Q ^ ^ £ < ^ ^ fl// ^ £ ^ 

2 f v,.M\ 1 and f(x) > c 2 for all x E C/B(o, ^r) where ci,c 2 > are 



- - . fa 2 /ii.M\l and fix) > t 

P(F-E[F]> U )<cxp|- — (15) con i^- n 



for M > and r * m\ I 1 1 1 1 



P(F-E[F]>a)<exp<^-- T ^ (16) 

I 2u J and F^(f,T) = Q(T),T ^ 00 for n> 2. 

for M = w/iere = (1 + i) ln(l + i) - i. Proo/; See |8] (sec tion 7.4, lemma 45). 



Lemmas 



Lemma 24. Denote Vd = r ^d +1 ^ the volume of a ball of 
Lemma 19. w a Poisson point process of intensity mdius 1 in R d flwt/ y J = d y d 2 f/zen 

measure ni\ dx for all t. 

R k+d 



Proof: Consider the point process Y^j,>i^Xi+Mi(t)- By / I 2 '! ^ x ~ ^ d f~ — j 

the displacement theorem ( flT4l ). it is a Poisson point process B (°^) 

of intensity measure dA t (i) characterized by: for all real k > —d. 

A t (A) = A f P(x + M(t) E A) dx Pr °° f: We haVe 

JR d 



= X dx p M (t)(y)l-{x+yeA} dy Jb(o.r) 



R 

x\ k dx = I r k d (V d r k ) 



= A / p M (t)(y) dy / l {x+y€A} dx 

JR d JR d 

= \k(A). 

Thus, it is a Poisson point process of intensity A dx. Now 
by thinning property, ujf 1 = T li >i 1 {h(t)=i}^X i +M i (t) is a 
Poisson point process of intensity A d x. ■ 



= K 



V d r k+d ~ 1 dr 
iR k+d 



d k + d 



